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Abstract 

We present a natural scenario for obtaining an open universe {Qq < 1) through 
Y^- inflation. In this scenario, there are two epochs of inflationary expansion — an epoch 



Q I of 'old inflation,' during which the inflaton field is stuck in a false vacuum, followed 



by an epoch of 'new inflation,' during which the inflaton field slowly rolls toward its 
true minimum. During the first epoch, inflation solves the smoothness and horizon 
problems. Then an open universe (with negative spatial curvature) is created by the 
nucleation of a single bubble. In effect Q is instantaneously 'reset' to zero. During the 
subsequent 'new' inflation Q rises toward unity. The value of Q today is calculable 
in terms of the parameters of the potential, and we show that obtaining values sig- 
nificantly different from zero or unity (though within the range < i7 < 1) does not 
require significant fine tuning. We compute the spectrum of density perturbations by 
evolving the Bunch-Davies vacuum modes across the bubble wall into its interior. 

[October 1994 — Revised: February 1995] 



1. Introduction 

According to conventional lore, inflation predicts a spatially flat universe with 
the cosmological density parameter Qq being very close to unity. It is therefore inter- 
esting to inquire whether inflation would be ruled out as a viable cosmological theory, 
or at least rendered sufficiently contrived to lose much of its underlying beauty and 
attractiveness, if it were to become observationally established that we live in an open 
universe with r2o < 1 today. In this paper we present a class of inflationary models 
that generate such an open universe in a natural way. This means that observations 
indicating Qq is less than unity can only limit the class of possible inflationary models 
and cannot rule out the whole idea. Our calculations shall show that in a sense the 
inflationary scenario is compatible with any value of Qq between zero and unity. 

Before presenting our scenario we flrst review the standard non-inflationary ar- 
gument against VIq ^ 1. For an equation of state of the form p ~ a~'^, (fi^ — 1) 
scales as a"'^'^, where a is the scale factor of the universe. This means that during a 
radiation dominated or matter dominated epoch, Q flows away from one. However, 
during an inflationary epoch (deflned here by the condition 7 < 2) i7 flows toward 
one. In the usual computation of the extreme unnaturalness of Qq being very different 
from unity, one evolves Q back to some much earlier time assuming that 7 ~ 3-4 and 
notes just how precisely Q must coincide with unity at that early epoch in order to 
evolve into a value of the order of unity today. If Qq is not exactly one, the deviation 
of this earlier value of Q from one gives a very small number, which is taken to reflect 
the degree of flne tuning required to obtain a nonflat universe. 

Such a calculation, however, is inappropriate if there is a preceding inflationary 
epoch, during which co-moving scales are pushed out of the Hubble radius, because in 
such a scenario the co-moving scale Xcurv characterizing the curvature of an open (or 
closed) universe is within the Hubble radius sometime during inflation, or perhaps at 
the beginning of inflation. If curvature is introduced as an initial condition prior to 
inflation, inflation acts to mitigate some of the unnaturalness of fi 7^ 1, but at a high 
price. By tuning the length of inflation to be short, one can by imposing curved initial 



conditions arrange to have significant curvature today, but in this case one would also 
have to put in by hand smoothness on scales of order Xcurv, and thus much of the 
initial motivation for proposing inflation would be lost. 

The natural way to obtain an open universe from inflation is to introduce spatial 
curvature during inflation, so that whatever inhomogeneities may exist prior to infla- 
tion are still erased, but spatial curvature is regenerated. In this way, inflation still 
provides a solution to the smoothness and horizon problems, while producing a non- 
flat universe. In such a scenario little flne tuning is required, because it is the length 
of inflation after the event introducing the curvature that must be adjusted to obtain 
the proper Qq. Said another way, one does not tune Qq but rather the logarithm of 
its difference from one, which requires little flne tuning. 

The remarkable fact that one can introduce spatial curvature within an inflating 
universe in a completely causal way, and without spoiling its homogeneity, was flrst 
revealed in the calculations of Coleman and De Luccia , prior to the invention of 
inflation. Their calculations showed that the spacetime interior of a (flnite) nucleating 
bubble actually encompasses an (inflnite) open FRW universe. The idea that a real- 
istic open universe could emerge from de Sitter space was flrst proposed by Gott in 
1982, prior to the introduction of the new inflationary scenario, and independently 
of Coleman and De Luccia's work. 

A realistic open universe scenario may be obtained via this mechanism, by com- 
bining old inflation with new inflation. The universe begins in a false vacuum in what 



f We emphasize here that what we mean by new inflation is simply inflation that proceeds 
whilst a scalar field is rolling slowly down a potential (i.e. 'slow-roll' inflation). In its original 
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incarnation, the new inflationary scenario ' assumed the universe began in a state close to 
thermal equilibrium before inflation began — this was presumed to explain how the inflaton field 
became localized around a potential maximum. We shall make the much milder assumption 
that somehow a scalar field became trapped in a false vacuum over a sufficiently large region 
for inflation to begin. This could have happened as a result of 'random' initial conditions, in 
which case our scenario would be more akin to chaotic inflation, 'special' initial conditions, 
such as the Hartle-Hawking-Vilenkin proposals, or even through a state of thermal equilibrium. 
Once the scalar field becomes trapped, and the epoch of 'old' inflation begins, the details of 
how it became trapped are of course very quickly erased. 



might be described as the 'old inflationary epoch,' during which any pre-existing in- 
homogeneities are redshifted away. The smoothness and horizon problems are solved 
during this first epoch of old inflation. Then a single bubble nucleates, inside of which 
the entire presently observable universe is formed. After bubble materialization, the 
surfaces on which the infiaton field 0, which we shall for the most part assume to be a 
single real scalar field, is constant are surfaces of constant negative spatial curvature. 
[If one ignores gravitational effects and thus considers bubble nucleation in Minkowski 
space, is constant on surfaces of constant s^ = (i? — x^ —y'^ — z'^), which for s^ > 
are surfaces of constant negative spatial curvature.] It is these surfaces that will be- 
come surfaces of constant cosmic density, or constant cosmic 'time' in the usual sense. 
However, if one would proceed directly to an FRW post-inflationary epoch by reheat- 
ing near the bubble wall one would obtain Q somewhat less than one at re-heating, 
and the present value Qq would be extremely close to zero. Such a scenario is clearly 
unacceptable. Instead we propose passing to an epoch of new inflation after bubble 
nucleation, tunneling from the false vacuum onto a slow-rolling potential. To obtain 
an open universe one wants the Hubble radius shortly after bubble nucleation to be 
of approximately the same co-moving size as the present Hubble radius, so that Qq is 
near one. This involves adjusting the length of the epoch of new inflation. In models 
of new inflation designed to explain fig = 1 the length of new inflation is arbitrary, 
provided that it is long enough to erase the initial inhomogeneities. What happens 
at the end of inflation is observationally accessible, but the earlier part of inflation is 
not, because what happened within the Hubble radius at earlier times now lies hidden 
on co-moving scales much larger than the Hubble radius of our presently observable 
universe. In the inflationary scenario that we consider, the entire new inflationary 
epoch is observationally accessible, and part of the old inflationary epoch as well. 

We emphasize that none of the ideas we have discussed so far are in any way 
new. Indeed they follow precisely the historical lines of development of 'old' and 'new' 
inflation. However we believe that the predictive nature of the 'single-bubble' scenario 
in a situation where the spatial curvature inside the bubble is not completely wiped 
out by 'new' inflation inside the bubble has been overlooked. It is this signiflcant 



lacuna in the literature which this paper seeks to fill. 

The bulk of this paper deals with calculating the spectrum of density perturba- 
tions in the open infiationary model just described. On 'small' co-moving scales — 
where small is defined relative to the scale of spatial curvature — the spectrum of 
density perturbations is approximately scale invariant, just as in conventional 'new' 
infiation. This is because small scales leave the Hubble radius at late times dur- 
ing the epoch of new infiation, when the curvature scale is much larger than the 
Hubble radius. Hence curvature is of little relevance in determining the small-scale 
perturbations. 

By contrast, on larger co- moving scales scales negative spatial curvature becomes 
increasingly relevant, and we expect deviations from scale invariance to occur. For 
large scales quantum fiuctuations that arise during the epoch of old infiation are also 
relevant. There are also other effects from the change in d'^V /dcjP' on passing through 
the bubble wall that lead to parametric excitations. These may be called 'moving 
mirror' effects ' . 

The zeroth order classical expanding bubble solution is 5*0(3, 1) invariant and de- 
scribed by a coordinate system that divides spacetime into three regions, as indicated 
in Figure 1. In region I the metric has the form 

ds'^ = -dt^ + a^{t) ■ [di^ + sinh2[^] rffi2^J (1.1) 

where dVt^^^-, = dO'^ + sin OdcjP'. The metric (1.1) represents an open expanding FRW 
universe. An important feature of this metric is that it is nonsingular - the usual 
'big-bang' singularity at a{t) = is here just a coordinate singularity - at small times 
a{t) ~ t and the metric [(1.1)] describes the Milne universe (see e.g. ref. 5), just a 
coordinate rewriting of Minkowski space. In region II the metric has the form 

ds'^ = da^ + 62(a) ■ [-dr^ + cosh2[r] rffi^g)] (1-2) 

and b{a)/a ^ 1 as o" — i> 0. In region HI, which essentially is an open FRW universe 
expanding backward in time, the metric the same form as in eqn. (1.1). In the perfect 



expanding bubble solution the scalar field is constant on surfaces of constant t and 
constant a as well as on the light cone of the origin, which separates the three regions. 

Only the upper half of the classical solution, sketched in in Figure 1, is physically 
realized. The lower half, which represents a bubble expanding backward in time, 
should be replaced with an instanton, known as the Coleman-de Luccia instanton, 
representing the classically forbidden process of bubble nucleation, as sketched in 
Figure 2. The instanton and the classical expanding bubble solutions are trivially 
related to each other by analytic continuation. 

At first it may seem paradoxical that an infinite open universe can fit inside an 
expanding bubble, which at any moment would seem to have a increasing finite (but 
never infinite) size, increasing at approximately the speed of light. Indeed if one 
were to use a time foliation in figure 2 in which horizontal lines would correspond to 
constant time hypersurfaces, at any finite time the bubble would have a finite size, 
and an observer using such a coordinate system would be quite justified in saying that 
the size of the bubble is always finite. But such a foliation would not correspond to 
what we today would consider a natural definition of constant time. For us constant 
cosmic time corresponds to a constant CMB temperature. To us a foliation that 
makes the spatial homogeneity and isotropy of the universe manifest seems preferred, 
and such a time slicing corresponds to the hyperboloids skteched in region I of figure 
2, which are infinite. This is how an infinite open universe can fit inside a bubble. 
Despite their contradictory appearance, the two points of view are quite compatible 
with each other. 

Another possible source of density perturbations, which we do not consider here, 
arises from imperfections in the initial bubble. In the semiclassical limit h ^ 
the scalar field configuration that results from quantum tunneling is fixed, with no 
random variations from the perfect expanding bubble solution. But as one considers 
corrections of higher order in h, small random variations occur, because the field does 
not tunnel exactly along the configuration of least Euclidean action, but rather along 
a nearby path. We believe that such perturbations in the tunneling process itself are 



insignificant, for two reasons. First, if the bubble nucleation is a strongly suppressed 
process, the fluctuations around it are likely to be even more strongly suppressed. 
And second, because in the subsequent classical expansion of the bubble the details 
of how the bubble was initially formed are likely to be quickly erased. The dominant 
sources of perturbations rather are the quantum fluctuations of the fields in de Sitter 
space prior to bubble nucleation, which pass across the bubble wall into its interior, 
and those that evolve later inside the bubble. 

Open inflation has been considered more recently in refs. 7-9. Lyth and Stewart, 
and Ratra and Peebles ' have considered density perturbations in an open inflation- 
ary scenario, a scenario of new inflation with negative spatial curvature put in by 
hand. They do not base their calculation of the spectrum of density perturbations 
on a speciflc model for the creation of negative curvature at the beginning of their 
inflationary scenario. Rather the discussion is conflned to what we call region I, and 
as initial conditions for the fluctuations of the inflaton fleld at t = 0, they impose 
the requirement that operators corresponding to the modes with asymptotic behavior 
t^*''^^ (as opposed to t'^^'^^^) as t ^ 0+ annihilate the vacuum. This state is some- 
times referred to as the 'conformal vacuum', and it has the following unappealing 
property. The energy density in scalar fleld fluctuations in their assumed quantum 
state actually diverges as one approaches t -^ 0+, which means that the coordinate 
singularity in the background solution at t = becomes a real singularity in this state. 
Thus one loses one of the principal attractions of open inflation, that the standard big 
bang singularity is removed. In contrast, with our choice of initial state, the energy 
density is nonsingular as t ^ 0+, and continues smoothly through the coordinate sin- 
gularity at t = 0, into de Sitter space. So whereas Lyth and Stewart, and Ratra and 
Peebles, have, for consistency, to simply assume their universe begins on a spatially 
homogeneous and isotropic slice some short time after t = 0, and thus give up on 
the inflationary solution of the horizon/homogeneity puzzles, we have as our starting 
point a long era of 'old' inflation outside the single-bubble universe, which in which 
those problems are solved. 

The quantum state of a scalar fleld of constant mass outside an expanding bubble 



in Minkowski space was considered by Rubakov, , by Vachaspati and Vilenkin, and 
Garriga and Vilenkin , using the thin-wall approximation. Sasaki et al. considered 
the quantum state inside such a bubble in the usual thin-wall approximation. Recently 
they have extended this work to the case of de Sitter space, but where the interior of 
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the bubble is empty, so no 'new' inflation and reheating occurs. 

The organization of the paper is the following. In section II we give a simple 
discussion of how the value of Qq depends on the parameters of the inflaton potential 
V{(f)) and the reheating temperature. In section III we give a brief review of bubble 
nucleation in de Sitter space in the semiclassical limit, in which there are no pertur- 
bations. The Euclidean bounce is the zeroth order solution, about which we consider 
small perturbations treated as free quantum fields. Section IV contains a new 'thin 
wall' approximate bounce solution which we use to show that there is a well defined 
region of parameter space in which significant simplifications of the nucleating and 
expanding bubble solutions occur. The analytic continuation of the Euclidean bounce 
into a spacetime with a Lorentzian spacetime signature possesses an 5*0(3, 1) sym- 
metry which we exploit to simply our calculations. Section V develops the formalism 
necessary to do this. In section VI we consider initial conditions prior to bubble nu- 
cleation. This involves rewriting the two-point function for the Bunch-Davies vacuum 
in terms of the region II mode functions. In section VII we deal with the evolution 
of cosmological perturbations in region I. In section VIII we deal with re-heating 
and observational consequences, and finally in section IX we present some concluding 
remarks. 



2. Calculating Q from Inflation 

The scenario we shall explore in this paper is based on a scalar potential V{(f)) 
of the general form illustrated Figure 3. We assume that the scalar field became 
stuck in a 'false vacuum' F in a region of the universe large enough for an extended 
period of infiation to occur. Once this has happened, it is very plausible that the 
universe is to a very good approximation described as de Sitter spacetime. This state 
of affairs is interrupted by the nucleation of a bubble. In a localized region of space 
the inflaton field tunnels quantum mechanically through the potential barrier from 
its initial value (pi to a value 0n, the value of the field at the center of the nucleated 
bubble at the instant of nucleation. We describe the nucleation process semiclassically, 
using a Euclidean 'bounce' solution. Inside the bubble, the field rolls down the 
potential toward (pr, where oscillations of the scalar field reheat the universe and lead 
to a conventional radiation dominated universe A spacetime picture of the bubble 
nucleation and subsequent reheating is shown in Figure 2. = 0„ on the future light 
cone of the point at the center of the nucleated bubble. Surfaces of constant within 
this light cone are spaces of constant negative curvature (open universes). They are 
also surfaces of constant density in the subsequent evolution of the universe, which 
we would call surfaces of constant cosmic time today. But these surfaces are of course 
infinite. So inside a finite bubble one has produced an infinite open universe. This 
remarkable picture was first revealed in the calculations of Coleman and de Luccia. 

We shall begin with a simple point. In this scenario the initial conditions prior 
to bubble nucleation and 'new' infiation are very definite — the universe has become 
very homogeneous and the value of has been fixed. From this follows a great deal 
of predictive power. In particular the value of Q today is not a free parameter; it is 
fully determined by the form of the potential, and by the microphysical details of the 
reheating process. So one can sensibly ask the question of how much 'fine tuning' 
of the parameters in the potential is needed in order to produce a given value of Q 
today. This is easily calculable with some reasonable simplifying approximations. We 



f The detailed structure of V{(j)) around 1/)^ shall play little role in our discussion. 



shall assume that the potential well around 0j is narrow, so that (f)n ~ (pi. We further 
assume that for most of the region between 0„ and 0^, the linear approximation 

V(0) ^ -/i30 (2.1) 

is valid. (This is easily generalized to an arbitrary power law form for V^(0), with 
very minor changes in the conclusions). For simplicity we have defined the value of 
at reheating, (pr in Figure 3, to be zero. As in all inflationary scenarios, we have to 
adjust the potential so that V = Q (no cosmological constant) in the universe after 
reheating. The motion of is described by the equation 

dV 

^ + ^H^ = -^ (2-2) 

with H as usual the Hubble constant, and we take to be homogeneous on surfaces 
of constant cosmic time t. The Einstein equation in an open universe is just 

\aj 3 2 a^ 

where a{t) is the scale factor. These equations describe the classical field solution for 
the interior of a nucleating bubble (see Figure 2). 

In the 'slow-roll' approximation, which becomes very good after a short initial 
transient, the term in ((2.2)) and the 0^ term in ((2.3)) may be dropped. The scale 
factor a{t) is then well described by an adiabatic approximate solution to (2.3): 



a(t) ^ — sinh 



dt B 



,2 87rG\/(0) 



B^ = ^-^. (2.4) 



At early times, the spatial curvature dominates and a{t) oc t, but as inflation sets 
in a{t) starts growing exponentially, with a{t) oc eJ . Once this happens, the 



curvature term in eqn. (2.3) becomes negligible and H k, B. The number of e- 
foldings during inflation is approximately 



tt2 

dt H = I H^^ " ' '^^ 



VM 






(2.5) 



Thus the number e-foldings depends only on the total change in the scalar field 
during inflation, and not on the parameter determining the slope of the potential, 
which we call /i. (The magnitude of the density perturbations produced during infla- 
tion does depend on yU, and requiring them to be small requires that /i <^ fnpi). 

In this scenario the value of Q begins at zero on the null surface defined by t = 
(or a = 0) but rapidly approaches unity during the inflation inside the bubble. From 
eqn. (2.3) one has 

where the last expression holds at reheating, and follows from eqn. (2.5). After 
reheating, we have il~^ — 1 oc (pa^)^^, growing as a^ in the radiation and as a in the 
matter eras, respectively. The value of Q today is then given by 

^0 ^ Y7^— 8^GA^5 ^^(^^j^ (2.7) 

where Tq is the GMB temperature today, T^-h that after reheating and T^q that at 
equal density of matter and radiation. If T^ is of order the electroweak scale, ^100 
GeV, then ^ ^ lO^^, but if Trh is of order the GUT scale, then A ^ 10^°. The 
qualitative behavior of eqn. (2.7) is easily understood. If A0 is small, then little 



t We have here assumed instantaneous reheating. If reheating is slower, there can be an in- 
termediate, matter dominated stage during which fl~^ — 1 scales inversely with temperature 
rather than temperature squared. This is a minor numerical detail. 
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inflation occurs and Qq is very close to zero. Conversely, if A0 is large, a lot of inflation 
occurs, and Qq is very close to unity. Generalizing the calculation to arbitrary power 
law form for the potential, V{(f)) oc 0", one finds that the factor Svr is replaced by 

8n/n. 

Defining the value of at reheating to be zero, the value of (pi (or equivalently of 
0n, because we have assumed they are nearly equal) required to obtain Qq today is 
given by 

A 



-^^ In 



fio'l - 1 



(2.^ 



»7r 

As a concrete numerical example, let us assume reheating up to the electroweak scale, 
so that A = 10 , and a linear potential, so that n = I. A value of r2o between .1 and 
.9 today then requires 1.48 < {(pi/rripi) < 1.54. 

Following this exercise, it seems to us that it is hard to argue that inflation is 
incompatible with an open universe. In any fundamental theory, one expects that the 
value of 0j would be fixed, presumably in Planck mass units. It would not seem very 
implausible that (pi should lie within the few percent of parameter space required to 
produce an interesting value of Qq significantly different from unity today. 

3. Background Solution 

The details of bubble nucleation in a curved space of nonvanishing constant cur- 
vature (i.e. de Sitter space) were first worked out by Coleman and de Luccia, 

[18l 

who generalized earlier work on bubble nucleation in flat Minkowski space, where 
gravitational back reaction from the nucleating bubble is not taken into account. 

Before outlining the modifications necessary to include the effects of gravity, we 
give a lightning review of bubble nucleation in Minkowski space through quantum 
tunneling at zero temperature. The bubble nucleation rate is 



A exp 



--Se [M^)] 



(3.1) 



where Se is the Euclidean action and (pti^) is the 5*0(4) symmetric Euclidean bounce 
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solution, satisfying the equation 

(920b 3 d(j)b _ dV 
ds^ s ds dd) 



(3.2) 



where s = (x^ + y"^ + z"^ + t^Y''^. The constant A has dimensions of [mass)'^^ and is 
harder to calculate, but may be estimated crudely as m^ where m is the mass of the 
scalar field in the false vacuum. 

The Euclidean bounce solution should be interpreted in the following manner. 
To materialize a bubble, the scalar field must pass through a classically forbidden 
region, under a classical potential barrier. For the Euclidean bounce given above 
(with its materialization center located at x = t = 0), the classically forbidden 
evolution occurs during t < 0, and is described using the Euclidean path integral 
formalism. The classical bounce represents the path of least Euclidean action to the 
classical turning point. For t > 0, the evolution is classically allowed, and thus treated 
classically. Formally one changes the signature of the metric from Euclidean for t < 
to Lorentzian for t > 0. In the Minkowski space case, one can regard this as distorting 
the path integral over (j) into the complex plane, to run over a (complex) saddle 
point. 

We now review the generalization with gravity taken into account. When 
^[0/1-] > 0) the underlying spacetime is de Sitter space, whose Euclideanized ver- 
sion is S*^ because under analytic continuation the metric ds^ = —dt^ + cosh [t] dVl?^o\ 
becomes ds^ = +dt^ + cos2[t] di^'jo\, which describes a four- dimensional sphere of 
radius R = H^^. 

The equations of motion for the Coleman-de Luccia bounce are derived by con- 
sidering the Euclidean action for the scalar field coupled to gravity, where 5*0(4) 
symmetry is imposed to restrict the form of the metric to 

ds'^ = da'^ + b'^{a) dnf^y (3.3) 
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They are 



^"(a) 



6V) 
b{a) 



Sb'ia) dV 



2 



^TtG f 1 ^, 2 



62((t) 3 12 



-0' ^ - \/[0(a)] 



(3.4) 



subject to the boundary conditions 



h{a = 0) = h{a = (Tmax) = 0, 

0'(o- = 0) = 0'(o- = (T^aa;) = 0. 



(3.5) 



(The latter two conditions follow from the continuity of the solution - the parameter 
a is related to the time t inside the bubble by cr = —it, so we have at sigma = t = 
that ((90/(9(7) = i{d(j)/dt), with 0(t) and 4>{cr) real.) 

In the flat space case (where the coupling to gravity is turned off), (Jmax = +00 
and the bounce action starts in the false vacuum, with 0((T = +00) = (j)jy. Because 
of the flrst derivative term (which behaves as a sort of friction term), V^[0(o" = 0)] is 
slightly higher than K[0j„]. With gravity turned on a^max is flnite, and the Euclidean 
instanton does not quite start in the false vacuum but rather partially into the barrier. 
As the coupling to gravity is further increased, the Coleman-de Luccia instanton 
approaches the Hawking-Moss instanton — in other words, (j){a = 0) and 0(o" = cr^ax) 
approach the maximum of the potential '' . For a discussion of the interpretation of 
the Hawking-Moss instanton, we refer the reader to refs. 2. 

If one continues the Coleman-de Luccia bounce to a spacetime with a Lorentzian 
spacetime signature, the functions 0((t) and b{a) describe the solution in region II. In 
region I (inside the bubble), the fleld rolls toward the true vacuum. In region IV the 
fleld rolls toward the false vacuum, undergoing decay and an inflnite number of oscil- 
lations. [In the maximal Lorentzian extension, regions III and V mirror the behavior 
in regions I and IV, but these regions lie in the domain of Euclidean evolution.] 

We shall be primarily concerned with the case where the size of the bubble is 
small so that a good deal of the Euclidean bounce is very nearly de Sitter space. In 
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this case we may use the Bunch-Davies vacuum (for a scalar field of constant mass in 
de Sitter space) as an initial condition for the fiuctuations in the scalar field at t = 0. 

There are two unsettling aspects of the Coleman-de Luccia instanton. The first 
is the global character of the tunneling process. All of de Sitter space is involved. 
There is no asymptotic region in which the deviations in the fields from their false 
vacuum values vanishes. Formally, this is the result of constructing the instanton from 
a closed geometry (i.e., choosing a time variable for analytic continuation such that 
the constant time hypersurfaces are compact). One would prefer an instanton based 
on the 'fiat' coordinatization of de Sitter space, in which the tunneling event would 
occur within a Hubble radius, and in which the instanton would have a tail falling 
off rapidly outside the Hubble radius. Such an instanton would seem more physical, 
but unfortunately no such instanton is known. The second unsettling aspect is that 
these instantons involve continuing the spacetime metric as well as the scalar field. 
In a sense one is considering complex metric configurations. So there is implicitly 
some sort of integration of spacetime metrics being performed (i.e., a path integral 
for quantum gravity), which one is then distorting to run over a complex saddle 
point. We shall in this paper adopt a conservative view, namely imagine that we are 
working in the regime where the nucleated bubble is much smaller than the Hubble 
radius H~^, so that the gravitational effects are small. Our bubble solutions are then 
well described by the Minkowski space solutions, and so (hopefully) little affected by 
these quantum gravitational issues. 

4. A New 'Thin Wall' Limit 

This open inflation scenario can be most simply explored in the following regime: 

(1.) The bubble nucleation rate is small, so that it is reasonable to suppose that 
the universe became very accurately described by de Sitter space, and the quantum 
fleld mode functions became well approximated by the Bunch-Davies vacuum modes. 

(2.) The size of the nucleated bubble is small compared to H~^. This allows us to 

14 



treat the nucleating bubble as a geometrical 'point' on the scales ~ iJ ^ of primary 
interest, and to suppose that it has very little effect on smaller scales. 

(3.) The bubble wall is thin. This allows us to match modes across a geometrical 
null surface, rather than evolve scalar field modes through the complicated structure 
of a bubble wall. 

(4.) The potential V{(f)) around the region 0„ is accurately approximated as linear. 
This again affords a useful technical simplification. Since we wish 'old' inflation to 
be followed by 'new' or 'slow roll' inflation it seems reasonable to assume that V{(f)) 
is a slowly varying (i.e., approximately linear) function around 0„. 

(5.) The vacuum energy in the false vacuum F is nearly equal to that along the 
light cone n. This requires that AV <^ V (see Figure 3). Again this allows us to 
treat the background de Sitter geometry as fixed in the matching of modes across the 
bubble wall. 

It is not obvious that such a regime exists, and it is the purpose of this section to 
show that it does. The 'thin wall' solution that we are looking for is not the standard 
one, described for example by Coleman and De Luccia, which holds when the false 
vacuum F is nearly degenerate with the true vacuum. Clearly we do not want to be 
in this regime, since then we would not have much 'new' inflation inside the bubble. 
In this section we find an analytic solution for the Euclidean 'bounce' and classical 
Minkowski bubble which holds in the regime we want. 

We consider the case where the mass m of the field in the false vacuum F is 
large so the potential well around 0j becomes very narrow, but its height AV" is fixed 
(see Figure 3). Our solution is exact in the limit as m becomes infinite. This is not 
in itself a case of much physical interest, but we do expect the infinite m solution 
to be a reasonable approximation to the case we are most interested in, where m is 
substantially larger than H. 

Anticipating that the final bubble size will be much smaller than H~^, we use 
eqn. (3.4), in the limit Ha ^ 1. The equation for b{a) tells us that at small a, 
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b r^ a r^ r, the usual radial variable in flat space. The scalar fleld equation then 
reduces to the flat spacetime equation for the Euclidean bounce: 

0"+V = -v:0 = V- (4.1) 

Here we have assumed the potential is linear over the entire relevant range. This 
is true right up to the beginning of the 'well' around (pi, which we assume to be 
very narrow. The bounce starts with = 0„ at r = 0. evolves according to eqn. 
(4.1) up to the edge of the 'well', and we have 

= 0n - ^r2. (4.2) 

As usual, the solution 0(r) may be thought of as the trajectory of a particle moving 
in a potential —V{(f)) with r playing the role of time. With this interpretation, the 
solution of the equation is then clear. accelerates according to (4.2), and hits the 
barrier —AV at some value of r we shall call r^,, its kinetic energy being converted 
into potential energy. (When the potential well is very narrow, is decelerated to 
a standstill in a very small interval of r, during which the loss of energy from the 
damping term is insigniflcant). Thus one flnds using eqn. (4.2) and setting = 0j 
that at r^, 

l0'2 = (0,-0,)^ = AK (4.3) 

This equation determines 0n, and then eqn. (4.2) determines r^ so that 



+ 4^-, n=J —. 4.4 



The Euclidean action for this solution is: 



Sj, = 2n' I r^dr\^-<P'^ + (V(0) - V(0.))] = ^^'^- (4-5) 







Let us check the consistency of this solution. The semiclassical approximation 
is only likely to be reliable in the case where the Euclidean action is large, which is 
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true if we take AV ~ /i^. The mass m has to be large in order that we can treat the 
false vacuum potential well as narrow — approximating it as (const.) + m^^cj) — (pr)'^, 



we require that the width of the well |0 — 0r| ~ ^/AV/m be much less than 0^ = 
4:AVfi~. Taking AV ~ /i , we see the narrow well approximation should hold well 
ii m ^ fi. Next we compare the bubble size r^ ~ yU"^ (Eqn. (4.4)) to the Hubble 

1 i 3 

radius H^^ ~ SrUpiV^^ ~ m^^n'^ ^ ^~i if /i ^ rupi. This justifies our use of the 
flat spacetime approximation in eqn. (4.1). 

These results mean that conditions 1 and 2 above are satisfied. Condition 3 
follows after noting that the classical solution for the expanding bubble is given by 
exactly the same solution (analytically continued) as the Euclidean bounce outside 
the forward light cone of the nucleated bubble's center, so conditions 3 and 2 are 
equivalent. Condition 4 is satisfied by construction, and condition 5 holds as long as 
AV <^ V, which is also true if yU ^ rUpi (recall from section 2 that the total range of 
from 0j to (f)r has to be of order rUpi to get a reasonable value for Qq). 

We have shown that there exists a regime in which the nucleating bubble, and 
the bubble wall become 'geometrical,' and the problem simplifies. This is the case if 
the mass in the false vacuum m obeys m^ ^ H^. However this case appears to be 
technically more complicated than the case where m^ = 2H^, for which a minimally 
coupled scalar field has the same dynamics as a conformally coupled, massless field. 
So for the remainder of the paper we shall restrict ourselves to m^ = 2H^, but 
nevertheless treat the bubble wall as 'thin'. We do this simply in order to be able to 
pursue the calculation right through to the end in a simple way. We do not anticipate 
any difficulties of principle in extending the treatment to m ^ H (see note added). 

The m^ = 2H^ case may actually be a reasonable approximation in the case of 
a two-field infiation theory, in which the first field, undergoes a bubble nucleation 
transition, in which it has a large mass squared in both phases, and there is no 
'slow roll' phase. If this couples to a second field, in such a way that the latter has 
m = 2H in the false vacuum, and m ^ after the tunneling event, then the second 
field can subsequently undergo a 'slow roll' transition. As far as the second field is 
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concerned, the bubble wall may then be treated as 'thin'. However we do not wish 
to enter into the complexities of two-field infiation here, and in fact regard it as a 
more interesting problem to extend the discussion to a single field, with m ^ H 
(see note added). 



5. SO {3, 1) Invariant Coordinates 

In this section we set up the coordinate systems and mode expansions most conve- 
nient for describing the expanding bubble and the perturbations about the perfect ex- 
panding bubble solution. The coordinate system that maximally exploits the 5*0(3, 1) 
invariance divides the (3 + l)-dimensional spacetime into three regions, indicated in 
Figure 2, in much the same way as Rindler coordinates divide (1 + l)-dimensional 
Minkowski space into four regions. Region I consists of the interior of the forward 
light cone of the origin O, region II consists of points with a spacelike separation 
relative to O, and region III consists of the interior of the backward light cone of O. 

In many parts of our calculations, we shall consider only the s-wave sector. Be- 
cause of the high underlying symmetry of the perturbations produced in the expand- 
ing bubble background, all of the density perturbations may be deduced from the 
s— wave sector. We calculate the two-point function choosing the origin as one of the 
two points, in which case only the s— wave sector contributes. Because of homogeneity 
and isotropy, the two-point function for all positions of the two points is determined, 
and thus the behavior of all other modes with higher angular momenta. 

For region I, the metric has the line element 

ds'^ = -dt^ + a^{t) [^^2 ^ gj^i^2 ^ ^^^2 ^ g.^2 g ^^2^ j ^ ^5_^^) 

and ^—g = a (t) sinh ^ sin 9. [For the special cases of de Sitter space and Minkowski 
space, a{t) = H^^ smh[Ht] and a(t) = t, respectively.] In region I, the equation of 
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motion for the scalar field of the unperturbed bubble solution is 



4>,{t) + 3^Mt) + v'[Mt)]=o (5.2) 

a[t) 



and the Einstein equation is 



*'*"- ' +^Uf^) ^VlUm^ (5^3) 



a2(t) a2(t) 3mli [2 \ dt 

For region II, the metric has the line element 

ds"^ = da^ + b'^{a) [-dr^ + cosh^ r {dO^ + sin^ 9 dcj)'^) ] , (5.4) 

and ^J—g = b^^a) cosh r sin 9. [For the special cases of de Sitter space and Minkowski 
space, b{a) = H^^ sin[iJo"] and h{a) = a, respectively.] Likewise, in region II the 
equation of motion for the unperturbed bubble solution is (denoting (9cr0 = 0') 

and the Einstein equation is 



6'2(or) 1 Stt I 1 fd(f)r^'^ 



6V)-6»'3,nJ,|2Va<r' -'''W-'lHO. (5.6) 

Note that these are exactly the same equations as those for the Euclidean 'bounce'. 
Also, eqns. (5.5) and (5.6) are related to eqns. (5.2) and (5.3) by the substitutions 
C, = T + iTT/2, t = ia, a = ib. 

The functions (poit) and (poicr) form the 'perfect' (i.e, 5*0(3, 1) invariant) bubble 
solution in regions I and II, respectively. If they are to match smoothly on the future 
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light cone of the nucleated bubble (a = t = 0), we must require that 



0o(^ = O) = 0o(t = O), 
0o(t = 0) = 0()(o- = 0). 



(5.7) 



We now consider, to linear order in 0, small perturbations of the perfect bubble 



+ 



(5.^ 



For region I, the d'Alembertian takes the form 



D 



-M^/^9'''d, 



-d^ + '-^d,. 



a{t) 
1 



d^ + 2 cothe dc + ^^ (dl + cot 6 00 + -k^dt 
sinh ^ V sin 9 



-d^ + 'Md-, 



a(t) ' a? it) 



9| + 2cothe9^ + — -2-(-L2) 
smh 4 



(5.9) 



Using separation of variables, we find the mode expansion for the solutions of the 
wave equation 

[n + m2(t)]0 = O (5.10) 

where m^(t) = V"[(j)Q{t)\. Writing 



0(t, e, d, 4>) = T{t- kh) i?(e; kh, I) YU9, 0), 



(5.11) 



one obtains 



a{t) 
(9| + 2coth^(95 



a2(t) 
/(/ + 1 



kl + m\t) 



sinh ^ 



kl 



T{t- kh) = 0, 
R{i-kh,l) = ^. 



(5.12) 
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For region II one has 



a^-al-mo. ' 



(9; + 2tanhr dr 



1 



b{a) "" b'^{a) [" ' ' ^' "' cosher' 

Writing 

0(a, r, e, 0) = S{a; uJas) Qi^; ^ds, l)YimiO, 0), 
one obtains the equations 

Sb'ia) „, r uj^ 



S" + 



s' + 



b{a 
Q + 2 tanh r Q 



d^ 
a) 



62 ( 

7(/ + l) 



m(a] 



cosh r 



2 ^^rfS 



5 = 0, 
Q = 0. 



Define 

and define 0(o", r; /, ?7i) similarly. 

In region I one may express the most general solution in the form 

oo 

^{t,^;l,m) = jdChR{^;Ch,l) 





X 



'4+\ChJ,m) T^+\t;Ch) + 4^\Ch,l,m) T^-\t;Ch) 



(5.13) 



(5.14) 



(5.15) 



(5.16) 



(5.17) 



We define ^^ = Cf + 1 for future convenience, and (+) and (— ) label two linearly 
independent solutions for the temporal evolution. Likewise, in region II we may 
expand the most general solution in the form 



+ CO 



(a,T:Lm) 



dCdS S{(t; Cds) 



X 



(a) 



A\yiCds;l,rn) Q^')iT;CdsJ)+A\j'iCds;l,m) Q^^^ir'XdsJ) 



Ml 



{5.U 



where ^^ = ^2 ^ ^^ 
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We now proceed to calculate the properties of the mode functions. Note that the 
hyperbolic spherical functions Ri{^; Ch) and the de Sitter functions Q^'^'{t] (dS^ ^^^ 
Q\°-){t] Qs, I) are universal, because they are completely determined by the 5*0(3, 1) 
symmetry; these functions do not depend on a(t), b{a), m?{t)^ and nn?{a). By contrast, 
the coefficient functions Ti{t;(h) and Si{(7;(ds) ^^^ ^ot universal and do depend on 
the functions a(t), b(a), m^(t), and m'^(a). 

Hyperbolic Spherical Functions. A nice discussion of the hyperbolic spherical 
functions with references to earlier work appears in ref. 17. These functions satisfy 
the equation 



^^ + 2cothe^ + (C^ + l)-A-^^^ 



Ri{^;Ch) = 



(5.19) 



where k? = Ch~^ -'-• -'■^ i^ convenient to rewrite (5.19) as 



sinh ^ 



d' 



d .,9 , /(/ + 1 



dcoshC,'^ ' ^ (icosh^ 

xRi{^;Ch) = 0. 



sinh ^ 



(5.20) 



It may be shown that the functions 



Ri{^;Ch)=Ni{Ch)-{-y^'-smh'^ 



HI . -^u/ 



rf^+1 



(icosh,^'+ 



1 cos(C/iO 



(5.21) 



satisfy (5.20), where 



NiiCh) 



¥M+^'m+^')---i(i+^') 



(5.22) 



and the orthogonality relation 



JdC sinh^e Ri{^;Ci) Rl{^;C2) = S{Ci - C2) 




(5.23) 
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and the completeness relation 



dCh Rli^a', Ch) Rli^b', Ch) = — . , o ^ — 

smli 4a 



(5.24) 



hold. Note that 



Roi^Xh) 



2 sin[ae] 
TT sinh ^ 



(5.25) 



de Sitter Functions. The de Sitter functions satisfy the equation 



d , d ,^ r 

--^ + 2tanhr — + (Cd5 



/(/ + !) 
cosh r 



Q/(r;Cd5) = 0, 



(5.26) 



which can be rewritten as 



cosh r 



rf2 



(i sinh r^ 



+ 3 sinh r 



d 



d sinh r 



+ (C 



dS 



/(/ + !) 
cosh r 



QKr;Cds) = (5.27) 



with c^^^c, = C(i5' + 1. We demonstrate that 



Qi {^; Cds) = miCds) cosh' r 



rf'+i 



d sinh r^+^ 



3±«CdST 



(5.2^ 



satisfies (5.26). 

Peculiar Mode Functions. We now discuss the other mode functions T;(t; (h) 
and Si{a;Qs), which depend on the particular choices for a{t), b(a), m^(t), and 
m^(a). In the general case, for the expanding bubble, the behavior of these func- 
tions depends on the particular choice of potential and will have to be evaluated 



t It is readily verified that Qo'''\T;Cds) = e'^'^'^ "■■'''/ cosh[T] satisfies (5.26) witli ^ = 0. We 
generalize this result inductively to all I. For this it is convenient to rescale z = sinhr and 
Qi{t) — cosh T wi{z), so that (5.26) becomes 



(z2 + l)w'i'{z) + {21 + i)zw[{z) + [(Cds' + 1) + l{l + 'i)]wi{z) ^ 0. 
If wi satisfies (5.29) for /, then w'l satisfies (5.29) for (/ + !). 



(5.29) 
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numerically. However, for certain special cases, we may solve for these functions 
analytically — in particular for a field of constant mass (i.e., whose mass does not vary 
with spacetime position) in de Sitter space and in Minkowski space. 

There are two reasons for considering these special cases. First of all, by consid- 
ering a massless field in Minkowski space, we may determine matching conditions at 
the light cone that are generally valid. More specifically, we would like to determine 
the matrix elements of the integral transform 



dCdS 




(5.30) 



T,^'''\Ch\Cds) Tj^^'^hChKds) "" [AfhCds;l,m) 



relating the expansion coefficients in regions I and II. 

This task would be quite simple if the individual mode functions in regions I and II 
did not become singular as t ^ and a ^ 0, respectively. One would simply require 
that and its transverse derivative be continuous across the light cone. However, 
each of the individual modes becomes highly oscillatory as one approaches the light 
cone. Therefore, wave packets that are well behaved near the light cone must be 
constructed as suitable smeared superpositions of individual modes. 

For matching the two expansions across the light cone, only the asymptotic be- 
havior of Ti{t;(h) and Si{a;((is) for small t and a is relevant. Therefore, we may 
calculate the matching for Minkowski space, with a{t) = t and b{a) = a, and the 
result obtained will be generally applicable, for all choices for the functions a{t), b{a), 
m^(t), and m'^(a), provided that the mode functions T/(t; (ii) and Si{a; (^g) are nor- 
malized to have the same asymptotic behavior near the light cone as the functions to 
calculate the matching in Minkowski space. 

Mode functions for Minkowski space. In Minkowski space, eqn. (5.12) be- 
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comes 



T + - T + 



m2 + ^ 



T = 0, 
which with the change of dependent variable T(t) = Y(t)/t becomes 



(5.31) 



Y + -Y 
t 



m 



t2 



F = 



(5.32) 



where /c^ = C| + 1- The solutions are Bessel functions of imaginary order J±ii^^(mt), 
and 



Tl^^\t■,Ch) = lJ±^cMt). 



(5.33) 



For the massless case m = 0, one has 

L 

Similarly in region II, for Minkowski space (5.15)b becomes 



(5.34) 



S" + -S' 
a 



UJ 



dS 



(7" 



m 



S = 0. 



(5.35) 



Making the substitution S{cr) = Z{a)/a, one obtains 



Z" + -Z' 
a 



CdS 2 

— ^ m 



Z = 



(5.36) 



where uj^g = CdS + 1- The solutions are modified Bessel functions of imaginary order, 
and 



S{a) = -K±i(^^^{ma). 



For small a 



S{a) ^ a 



±iC,dS-l - ^ 



±«Cdsln[cr] 



a 



(5.37) 



(5.3^ 



For rn^ > only one solution has acceptable behavior for large a. For the massless 
case, there are twice as many acceptable solutions, since both solutions oscillate as 
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a -^ oo. This doubling has a simple explanation. For the massive case, the light cone 
regarded as a Cauchy surface is equivalent to the Cauchy surface defined by r = 0, 
because for a massive field, at least massive in the limit a -^ oo, all of the information 
on the r = surface must eventually propagate into the light cone. For a massless 
field this is not so. There is a set of modes that propagate to null infinity ahead of 
the light cone, and hence never enter the light cone. 

Coordinate Systems for de Sitter Space. Before giving the mode functions 
for de Sitter space, to establish notation we write down the four sets of coordinates for 
de Sitter space used in the course of our calculations — which we shall call 'hyperbolic,' 
'flat,' 'closed,' and 'embedded' coordinates — and the transformations between them. 

One may construct (3 + l)-dimensional de Sitter space as an embedding in (4+1)- 
dimensional flat Minkowski space, defined by the equation 

x'^ + f + z'^ + v? -w'^ = 1. (5.39) 

We call the coordinates w, m, x, |/, z 'embedded' coordinates. 
The 'closed' coordinates are then given by 

w = sinh[trf], 

u = cosh[tc;] cos[xd], 

z = cosh[tc;] sm[xd] cos[9], (5.40) 

X = cosh[tc;] sm[xd] sin [6*] cos[0], 

y = cosh[tc/] sin[xd] sin [6*] sin[0]. 

They cover all of de Sitter space, and the line element is 

ds^ = -dtli + cosh^ [ta] ■ {dxli + sin2 Xdd^f^) ) • (5.41) 

In the 'flat' coordinates, the constant time hypersurfaces are null planes in the 
five-dimensional Minkowski space, into which the de Sitter space is embedded. More 
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specifically, 

tf = \n[w + u], 

f (5.42) 

w + u 
and the line element is 

rfs2 _ _^^2 ^ g2t^ . ^^^2 ^ rjdnf^^] . (5.43) 

where f = [x^ + y"^ + z^Y'"^ . The fiat coordinates cover only the half of de Sitter space 
defined by w + -u > 0. The other half is covered by another set of 'fiat' coordinates. 

Next we discuss the region II 'hyperbolic' coordinates 

w = sinh[r] sin[o"], 
u = cos[o'], 

X = cosh[r] sin[o"] cos[^], (5.44) 

y = cosh[r] sin [a] sin[6'] cos[0], 
z = cosh[r] sin[o"] sin[^] sin[0], 
in which the line element is 

rfs2 ^ ^^2 ^ gjj^2[^] . ^_^^2 ^ cosh2[r]rfl]^2)] ' (5-45) 

Finally we have the region I 'hyperbolic' coordinates 

w = sinh[t] cosh[^], 
u = cosh[t], 

X = sinh[t] sinh[^] sin[^] cos[0], (5.46) 

y = sinh[t] sinh[^] sin[^] sin[0], 
z = sinh[t] sinh[^] cos[^], 
in which the line element is 

rfs2 _ _^^2 ^ ginh2[t] . [rf^2 ^ sinh2[{]rffi2^J (5.47) 

These coordinates cover the forward light cone of the point P described by the em- 
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bedded coordinates u = 1, x = y = z = w = 0. To cover all of de Sitter space with 
hyperbolic coordinates, three more coordinates patches of a form similar to region I 
are required. Region III consists of the backward light cone of P. Region IV consists 
of the forward light cone of the antipodal point of P, which we call P, and which is 
described by the embedded coordinates u = —1, x = y = z = id = 0. Finally, region 
V is the backward light cone of P. 

Peculiar Functions for de Sitter Space. For de Sitter space in region II, the 
scale factor takes the form b{a) = H"^ sm[Ha]. For simplicity we set H = 1. Eqn. 
(5.15)a becomes 



S"{(t) + 3 cot [a] S'{a) - 77? S + 



sin \a\ 



S{a) = 0, 



(5.48) 



where < c < tt. Let x = cosa and S = F/sin[cr] = F/[l — x^)^'^, so that eqn. 
(5.48) becomes 



d 
dx 



X 



dF 
dx 



+ 



(2 - 771^) + 



Cds" 



X^ 



0, 



(5.49) 



which is the Legendre differential equation, with linearly independent solutions 



P 



±«CdS / 



^iJCdS 



(x) and QV^^^ix) with v{v + 1) = (2 - m^). 



For 771^ = 2, which we shall call the 'conformal mass' case (because the quadratic 
term in the action has the same effect as no mass with conformal coupling), we solve 
eqn. (5.49) by making the change of variable x = tanh[-u], with — oo < u < +oo, so 
that eqn. (5.49) becomes 



dv? 



{2 - 77?) sedi\u] F + Cds" F = Q, 



(5.50) 



which for m^ = 2 has the solutions F = e^*^'*^". For the massless case (m^ = 0) one 
has F = « - tanh[M])e±*^". 
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6. Initial Conditions 

One of the remarkable properties of a de Sitter space background is its ability 
to erase intial conditions through exponential expansion. After a sufficient number 
of e-foldings, the state of the scalar field within a Hubble volume becomes almost 
completely determined and almost completely independent of the intial conditions. 

In this section we discuss the quantum fiuctuations of the infiaton field prior to 
bubble nucleation, when the infiaton field is stuck in the false vacuum. There exist 
extensive discussions in the literature of the quantum fiuctuations in a free scalar field 
of fixed mass in de Sitter space. The principal result of this section is an expansion of 
the two-point function describing the Bunch-Davies vacuum (see ref. 5 and references 
therein) in terms of the region II hyperbolic mode functions. To follow the evolution 
of the quantum fiuctuations through the bubble wall, it is necessary to describe the 
initial state in this manner. 

The choice of the Bunch-Davies vacuum as a reasonable initial condition can be 
justified in the following way. In de Sitter space, if one follows any scalar field mode 
in the fiat spatial slicing, its proper wavelength begins inside the Hubble radius, and 
is exponentially stretched outside the Hubble radius as time proceeds. So the early 
evolution is effectively the same as in Minkowski spacetime. The energy density in the 
scalar field fiuctuations is the renormalized sum of the energy in each field mode - if 
this is to be finite, it is clear that the very short wavelength modes must begin in their 
ground state (i.e. the Minkowski vacuum). It follows that after a sufficient number 
of expansion times, the state of the scalar field (as observed over a fixed physical 
volume) becomes independent of the initial state of the scalar field, because all one 
sees are the modes which began in their ground state. This 'evolved Minkowski space 
vacuum' is the Bunch-Davies vacuum. 

For the two-point function in the Bunch-Davies vacuum \Obd) we shall use the 
Wightman function 

G^+\X,X') = (Ob^|0(X)0(XO|Ob^). (6.1) 
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where X = {w, u, x, y, z) are the 'embedded' coordinates of the previous section. Since 
G^'^> is invariant under the action of the connected part of 5*0(4, 1), we express G^^' 
in terms of the invariant 



/(X,X') 



-ww' + uu + XX + yy' + zz 



iee{X,X') 



(6.2) 



Note that /(X, X') > +1 when X hes inside either the forward or backward hght cone 
ofX'. /(X,X') = +1 when X hes on the hght cone of X'. When- 1 </(X,X') < +1, 
there exists a spacehke geodesic connecting X and X' . When J(X, X') = —1, then 
X hes on the hght cone of the antipodal point of X' . Finally, when /(X, X') < —1, 
X lies in the interior of the light cone of the antipodal point of X', and in this case 
there is no geodesic connecting X to X' . We add a small imaginary part to /(X, X') 
to indicate the time ordering of X and X' . e{X,X') vanishes when /(X, X') < 1, is 
equal to +1 when X lies inside the forward light cone of X', and is equal to —1 when 
X lies inside the backward light cone of X' . 

We have 



G^+\X,X') 



\ - v^\ 



IGtt 



• secivrz/l • F 



3 3 /X,X' +1 

- - z/, - + z/: 2: ^ '- 

2 '2 2 



(6.3) 



where z/ = | — vr? . Except for special degenerate values of z/, the hypergeometric 
function F[^ — z/, ^ + z/; 2; 2;] has a branch point at z = and a cut extends from 
2; = +1 along the real axis to 2; = +00. The infinitesimal imaginary part of / provides 
a compact notation for indicating on which side of the cut the hypergeometric function 
should be evaluated. We shall for the most part work in units where H^ = 1, and 
restore H by dimensional analysis later. 

To simplify the algebra, we shall in this paper only consider the special case where 



m 



2, for which the Wightman function has the particularly simple form 



G^+\X,X') 



87r2 /(X,X')-1' 



(6.4) 
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In this case the spatial mode functions S{cr; () also have a particularly simple form 

S(a:C) = r-r exp TiCtanh^ fcoso-)] = ; — -— (6.5) 

sm[o-] L ^ '-^ sech [u] 

where tanh[-u] = cos[o']. To normalize our mode functions, we use the nondegenerate 
bilinear form 



{f,g) = {-t)J rfS'^ ^f{X)[d,9{X)] - [d,f{XMX)j (6.6) 

E 

where S is a Cauchy surface, with unit normal n^, and rfS^ = dJ^n^, with dT, the 
volume element on S. For the s-wave, the normalized mode functions are 



-') 



e*Cw g=F«ICk 



4:TTy/( sech [u] cosh[r] ' 



which satisfy the relations 



(/^\4+^) = -5(C-0, 

(/r\4^^)=o. 



(6.8) 



We choose this normalization so that when we expand the s-wave component of the 
field operator as 

+00 

0(X) = [ dC \j<+\X; C)a(C) + I<-\X; C)at(C)l , (6.9) 



— oo 



the usual commutation relations [d(C),a^(CO] = '^(C ^ CO) [^(C))Q^(CO] = 0) ^^^ 
[a^(C),a+(CO] =Oresuh. 
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To expand G^^' in terms of the region II mode functions, we calculate the matrix 
elements 



M(Ci; ±i|C2; ±2) = f^^'\X; Ci) ox G^+\X, Y) oy f^^'\Y- C2) (6.10) 

where the contraction o is defined according to eqn. (6.6). We compute both products 
choosing the surface defined by r = to be the Cauchy surface S. 

In terms of the region II hyperbolic coordinates, 

/(X, X ) = cos[a] cosla ] 

+ sin [a] sinfo"'] (cosh[r] cosh[r'] cos[0] — sinh[r] sinh[r']) — ie sir — t') 
= tanh[-u] tanh[-u'] (6-11) 

+ sech [u] sech [u] (Z cosh [r] cosh [r'] — sinh[r] sinh[r']) 
— ie eij — T ) 

where Z = cos Q[9,(f);9', 0'] = cos 6 cos 6' + sin 9 sin 9' cos(0 — 0') . 
Consequently, 



g(+)(x,x'^ ^ 



Stt 



2 



X 



tanh[-u] tanh[-u'] + sech [u] sech [«'] (Z cosh [r] cosh [r'] — sinh[r] sinh[r']) 

-1 

(6.12) 

-1 1 



1 — ie sir — T ) 



Stt^ sech [u] sech [u'] cosh[r] cosh[r'] 

1 



Z — cosh[-u — u'] sech [r] sech [r'] — tanh[r] tanh[r'] — ie e{T — t') 
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We evaluate 



gW(X',X)o;,/(±)(X;C) 

+00 +1 



,-« 



du sech ^[u] {2n) / dZ 



47r2 sech [u] sech [u'] cosh[r] cosh[r'] 



(6.13) 



Z — cosh[-u — u'] sech [r] sech [t'] — tanh[r] tanh[r'] — ie s{t') 



sech [u] ^ ) I 47]-iy[^ sech [u] cosh[r] 



at r = 0. (Note primed and unprimed indices have been interchanged.) Define 



-1 +00 
H{u\t't) = [ dZ [ due+'^'' 



X 



Z — cosh[-u — u'] sech [r] sech [t'] — tanh[r] tanh[r'] — ie s{t') 

+00 



du e+*^"ln 



—00 

+CXD 



du e+*^"ln 



cosh[-u — u'] sech [r] sech [r'] + tanh[r] tanh[r'] — 1 + ie £{t') 
cosh[M — u'] sech [r] sech [t'] + tanh[r] tanh[r'] + 1 + ie sir') 

cosh[-u — u'] — cosh[r' — r] + ie e(r') 
cosh[-u — u'] + cosh[r' + t'] + ie ^(t') 



(6.14) 



so that 



gW(X',X)o;,/(±)(X;C) 
+i 



167r2^^ sech [m'] cosh[r'] 



-dH{u\ t', r) 



TmHiu',T',r) 



(6.15) 



at r = 0. 
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We now evaluate 



+00 



du e+*^" In 



+CX) 

=+^^"' / du e+^^" In 



cosh[-u — u'] — cosh[r' — t] + ie e{T') 
cosh[-u — u'] + cosh[r' + r] + ie e{r') 

cosh[-u] — cosh[r' — t] + ie s{t') 



cosh[-u'] + cosh[r' + t] + ie £{t') 



(6.16) 



Define 



r(±) 



du e+^^" In 



cosh[-u] — cosh[r' — r] ±ie 
cosh[M] + cosh[r' + r] ± ie 



(6.17) 



We first evaluate Jg . Examining the properties of the integrand in the complex 
plane, we find two branch points near the real axis, just above or below the two 
points u = ±{t' — t). The branch point to the left at -u = — |r' — r| +ie lies just above 
the real axis and the branch point at m = +|r' — t\ — ie lies just below the real axis. 
[Under the transformation +ie — > —ie this situation is reversed.] 

We now consider the contour C from — T to +T to T + 27ii to —T + 27ri and 
finally back to — T, considered in the limit T -^ 00. Inside the rectangle enclosed by 
the contour, in addition to the branch point at -u = —\t' — t\ +ie, there are three more 
branch points at M = — |r'+r|+i7r— ie, -u = +|r'+r|+i7r+ie, and-u = +|r'— r|+i27r— ie. 
One branch cut connects u = —\t' — t\ + ie to u = —\r' + r\ + in — ie, and another 
branch cut connects u = +\t' + r\ + in + ie to u = +\t' — t\ + i2n — ie. In the limit 
T ^ 00 the contributions from the vertical parts of the contours vanish. Since a 
translation by 27ri has the effect of multiplying the integrand by e^^''^. 
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du e+*^" In 



C 



cosh[-u] — cosh[r' — r] + ie 
cosh[-u] + cosh[r' + r] + ie 



— |t +r|+i7r 

:-27vi) 



I *»"^-"-p") / 



+ |r'-r|+2i7r 

m] I du e+*''" 

+ |r'+r|+i7r 



(6.18) 



27r 



Cg-iC|r'-r| _ g-C^r-iCk'+^l) _ (^-CT^+iClr'+r] _ g-2C7r+iC|r'-rl 



where the integral is evaluated by deforming the contour into two pieces surrounding 
the two branch cuts. It follows that 



Similarly, 



/(+) 



27r 



Csinh[C7r] 
27r 



cos[C|'r' + t\] — cosh[^7r — i(\T' — t\] 
(cos[C|r' + t\] - cosh[C7r]) cos[C|r' - t\] 



Csinh[C7r; 

+ i sinh[^7r] sin[^|r' — t\] 



(6.19) 



/(-) 



2tt 



Csinh[C7r] 
27r 



cos[C|'r' + t\] — cosh[^7r + i(\T' — t\] 
(cos[C|t' + t\] — cosh[^7r]) cos[^|r' — r 



Csinh[Cvr 

— isinh[^7r] sin[^|r' — t\] 



(6.20) 



It follows for both signs of r' that 



H 



27r 



,+iCu' 



Csinh[C7r] 
+ isinh[^7r] sin[^(r' — r)] >. 



cos[C(t' + t)] — cosh[^7r]) cos[C(t' — r)] 



(6.21) 
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Therefore 

(gWo/W)(X') 



o+iC,u' 



^7riy[^sinh[^7r] sech [m'] cosh[r'] 



X [cos[C('r' + r)] — coshfCvr] cos[C(t' — r)] + i sinh[Cvr] sin[C(r' — r)]] 



(6.22) 



Since 



j cosICt'I = Tii^ e 



eqn. (6.22) may be rewritten (setting r = 0) as 



tACW 



(g(+)o/(±))(xO 



1 



sinh[(^7r] 
1 



c 



-e^^l^l^' sinh[C7r] T ^ ■ {-e±''l^l^' + e^^l"^!^' cosh[C7r]} 



+iC,u' pT*ICk' 



4^v4Cl ^^^^ [""'] cosh[r'] 



X 



Therefore 



2sinh[|C|7r] 



G(+)o f ^^, 



T 



(+) 



e+ICkg^^l^l^' ± e^^'"^'^' 



/, 



(+) 



/. 



(-) 



-PW)-f '^. 



(-) 



-1 



=+ICk 



-1 



/, 



+ 



C 



e+|Ck_e-ICI^ 



.1 -e-\^\- ) \f^- 
and the normahzed positive and negative frequency mode functions are 



g 



(^^ _ elCk/2/(+) - e-ICIVV(-) 



(6.23) 



(6.24) 



(6.25) 



(6.26) 
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From the fact that G'-^-' obeys the wave equation everywhere in both variables and 
from the orthonormahty relations for the mode functions given in eqn. (6.9), it follows 
that 



X (±l) (T2) l^ft^Y) OyG^^\YX)oy, /f^)(rO) 



x/f^Hx)/f^)(X') 



The matrix elements /} (V) oy G^^'{Y^Y')oy' /) ^(V') are given in eqn. (6.25). 



7. Cosmological Perturbations 

In the previous sections, we calculated the quantum fluctuations of the infiaton 
field 0(a;) about the classical background solution 0b (x) treating the curved space 
as fixed background, determined by the background solution with no fiuctuations. 
To ignore the back reaction of the linear perturbations on the background spacetime 
is a valid approximation when the perturbations are well within the Hubble radius. 
However, as the expansion of the universe causes the size of a perturbation to increase 
relative to the Hubble length, this approximation breaks down. In this section we 
consider the coupling to gravity of the scalar field perturbations — in other words, how 
the quantum fiuctuations in the infiaton field translate into density perturbations as 
they are pushed outside the Hubble radius. It is permissible to treat the fiuctuations 
in the scalar field (and their associated metric perturbations) as classical fields. 



There are many physically equivalent ways to treat density perturbations. Here 
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we use the gauge-invariant formalism of Bardeen, closely following the notation in 
the review in ref. 31. We consider only scalar perturbations and write the metric in 
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the form 



ds = a {rf) 



-(1 + 2^)dTf + 2B\i dr] dx' + {(1 - 2ilj)-fij + 2E\ij]dx^ dx^ 



(7.1^ 



where rj is conformal time, and the spacetime functions ip, ip, B and E are hnear metric 
perturbations. Primes denote derivatives with respect to conformal time and Ti = 
a' /a = aH. Covariant spatial differentiation with respect to the spatial background 
metric •jij is denoted DjF = Fy etc. 



We set 



7y dx' dx^ = di^ + fiOdnf^) (7-2) 



where 




/(O = U, for /C = 0, (7.3) 



Although we are primarily interested in the open geometry (/C = —1), we consider 
the cases /C = —1, 0, and +1, (corresponding to open, flat, and closed spatial geome- 
tries, respectively) to allow comparison between the open and fiat cases, so that the 
influence of spatial curvature is manifest. In comparing open and fiat geometries, it 
should be k;ept in mind that spatial curvature affects the eigenvalues of the Laplacian 
V^ = 'j^'' DiDj used in the mode decomposition. For /C = 0, the eigenvalues of — V^ 
are /c^ where k ranges from to +oo; for /C = —1, the eigenvalues of — V^ are C^ + 1 
where ( ranges from to +oo. 



f In ref. 31 the metric is written as {dr^ + r^dQ%.)/{l + /Cr^/4)^. The substitutions r = 
2tanh[^/2] for /C = — 1 and r = 2tan[^/2] for /C = +1 demonstrate equivalence to the line 
element in cqn. (7.2). Likewise, in ref. 30 the line element is written as df^/(l — /Cf^)+f^dri?2)- 
Equivalence to eqn. (7.2) may be demonstrated by the substitutions f = sinh[^] and f — sin[^] 
for /C = — 1 and /C = +1, respectively. 
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From these linear perturbations, we construct the 'gauge-invariant' variables 

<^ = ^ + - [{B-E') a]', 



^ = ^ [B-E'], 



(7.4) 



where 'gauge-invariant' means invariant to linear order under infinitesimal coordinate 
transformations. We construct the 'gauge- invariant' perturbations of the Einstein 
tensor 

Sg^ = ^Gg + [ ^Gj]]' {B - E') 



5Q\ 



5G', + [ 'GO - i 'Gl] (B - E') 
5g'^=SG] + [^G)]' {B~E'), 



(7.5) 



and similarly of the stress-energy tensor 

5T,^ = 5TlJ + ['^T^nB-E'), 

5%' = 5Tf + [^Tf-^-'T^^]{B-E% 

5r^ = ST] + [^T]]' {B-E'). 



(7.6) 



In terms of the gauge-invariant metric perturbations, 



(7.7) 



where D = ($ — ^). The gauge-invariant stress-energy from the infiaton field is 



5T^ = ^ 



6Tt 



sr; 



+cp',^^-<P[{54>)' + a'V^{54>) 



(7.8) 



where (50) = (50) + 0^, (-B — E') is the gauge-invariant variable for the infiaton field 
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perturbation. The linearized Einstein equation 6Q'^ = (SttG) 57T gives D = (so 
that \E' = $) and the equations 

3(/C - 7i;2)$ - 3n^' + V2<l> = (AnG) [-0'^ ^^ + 0^, (50)' + a^V^ (50)] , 

H$ + $' = (4716-) [0';, (50)] , (7.9) 

$" + 37^$' + {2n' + n^ -K:)^ = {A-kG) ■ [-0;, 2$ + 0'^ (50)' - a^V^ (50)] . 

Because of the background Einstein equations 



3r.,2 



— [TT + K]= (SttG) ■ (+pb) = (SttG) 



^0'.' + n06) 



■ [2n' + n^ + ]c] = (SttG) ■ (-P;,) = (sttG) 



2a 



^ 0'6' + V^(06) 



(7.10) 



it follows that (47rG')0l = (7i — H' + fC). Thus eqn. (7.9) may be rewritten as 



7^$ + $' = (47rG') • [0;, (50)] , 



(7.11) 



We now subtract the first equation from the third equation. Because of the 
background equation of motion for the scalar field 



0;; + 27^0^ + a2\/^(0;,) = O, 



(7.12) 



we can subtract [47i — 2(0"/0^)] times the second equation, obtaining 
$" + 2(7^-%)$'- V2$ + (271' - 2n-^ - 4/C J $ = 0. 



(7.13) 



We rewrite eqn. (7.13) in terms of proper time, using the relation dt = a{r])dr], so 
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that 



$+ij_2^U + 4 (-V2 - 4x:) $ + 2 ij - //^ U = 0. 



(7.14) 



We now relate the gauge invariant potential perturbations to the gauge invariant 
scalar field perturbations. From eqns. (7.11) and (7.13) it follows that 






1 



1 



{A-kG) <P', 
-1 1 



- [$' + n^] , 



(AttG) 






C^ + 5 + n' - n^ ] $ 



(7.15) 



so that in terms of proper time derivatives 



dt{S4>) 



1 



-1 1 



(4vrG) 0, 






^+0.5,f^iK 



(7.16) 



To calculate the evolution of the density perturbations during the early part 
of region I infiation (until Q is close to one), we use a linear approximation to the 
potential. In the general case, at t = (at the end of the Coleman-de Luccia bounce), 
0b = and the potential slopes downward, toward the true vacuum. To make the 
problem tractable analytically, we take 

\/(0) = Vo- V^<P, 

where we set 0b(t = 0) = and take V(f, to be constant. We first solve for the evolution 
of (ph, for simplicity assuming that a{t) is well approximated by a{t) = sinh[t] until Q 
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is close to one. Thus the equation of motion for 0^ is 



3coth[t] <i)b = V^ 



(7.17) 



and the solution with 



at t = is 



Mt) = y,<t> 



cosh^[t] -3 cosh [t] + 2 
3sinh^[t] 



(7.1^ 



which for small t behaves as \V^(f,t and for large t as iV^^. It follows that 

hit) 3 



0j,(t) sinh[t] ■ (cosh[t] + 2) 



(7.19) 



so that eqn. (7.14) becomes 
1 



$ 



sinh[t] 



cosh[t] 



6 



cosh[t] + 2 



$ 



+ 



sinh^[t] 



C + 5 



4- (2cosh[t] + 1) 
cosh[t] + 2 



(7.20) 



$ = 0. 



To solve the eqn. (7.20) it is convenient to replace the variable t with the conformal 
time variable rj = ln[tanh[t/2]] , thus mapping the interval (0 < t < +oo) into 
(— oo < T] < 0), so that eqn. (7.20) becomes 



$' 



3 - e'^n 



iC' + 5) - 



4(3 + e^^) 



$ = 0. 



(7.21) 



:J 



The general solution to eqn. (7.21) is 



$ =c 



(+) 

+C(_) 



.+iCv . e'^ . ( 1 



-Kv . e^ . ( 1 



(C + 
3(C-0 

3(C + 



,277 



=277 



(7.22) 



For small t (i.e., 77 ^ —00), t ^ 2e'' and $ ~ C(-_|_-)e''e+*^'' + ct\e^e ^^'^ 



+ 



(t/2)(t/2)+*^ + C(_)(t/2)(t/2)-*^. For large t (i.e., 77 -^ 0-), $ is dominated by 



t We thank Bharat Ratra for pointing out that eqn. (7.21) can be solved analyticaUy. 
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the growing mode and 






*-^Wo 7^ +^(-)o(^^-)- (7-23) 



Finally, inserting the small t asymptotic forms $ ~ t^*^+\ a ^ t, H ^ t ^, 
(j)fj Ri jV^fjjt, and (ph ~ |V^0 into eqn. (7.16), we obtain the small t matching condition 

In spatially flat inflation the conserved quantity 

2 7^-1$' + $ 

[called ( in ref. 31 ] is very useful to track density perturbations on super-Hubble 
radius scales. During inflation 



w = pJpb = (-1) ■ -^-^ fS; ~ ^^^^ - 1. (7.26) 

Assuming the absence of entropy perturbations (i.e., that the perturbations obey the 
equation of state Pb{pb) given by the background solution) and neglecting the spatial 
derivative term in the evolution equation for the gauge-invariant potential $, one 
obtains x = 0. On super-Hubble radius scales these assumptions hold reasonably 
well. In flat inflation the existence of this conserved quantity provides an elegant 
way of demonstrating that the spectrum of density perturbations at late times is 
independent of the details of re-heating. The precise way in which w evolves, from 
a value slightly greater than —1 during inflation to g during radiation domination, 
and later to nearly zero during matter domination, does not affect the flnal density 
perturbations. 
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Unfortunately, the conservation of x for super-Hubble radius scales does not com- 
pletely generalize to an open expanding universe, x varies with time when Vt is not 
close to one, in our scenario during the early part of region I inflation and then much 
later, during the latter part of the matter domination, when the universe becomes 
curvature dominated again. Therefore we may use x i^i the following way. During 
the early part of inflation (when Vt is significantly less than one), we calculate the 
evolution of the modes explicitly. Then at a later time (before reheating but when Vt 
is still very nearly one), when all of the modes of interest are well outside the Hubble 
radius, we calculate x- (It should be noted that while Vt is less than one, x varies in 
time in a manner independent of wave number. Therefore, all modes are affected in 
the same way, and only the overall normalization of the power spectrum is altered. 
Its shape remains the same. This is because the evolution of density perturbations 
on super-Hubble radius scales is essentially a local process.) 

8. Power Spectrum of the Inflationary Open Universe 

In the last three sections we developed the tools for calculating density pertur- 
bations. In this section we put together the various pieces to give a concrete result. 
To keep the calculation as simple as possible, we consider the following idealized sce- 
nario. In region II we assign to the inflaton field perturbations a mass m^ = 2, which 
changes discontinuously across the light cone to m = in region I. This should be 
a reasonable approximation when the bubble wall is very thin and the bubble size is 
small compared to the Hubble length. [The scenarios with two scalar fields work best 
for this regime. See the discussion at the end of Section IV.] 

We start with the mode functions associated with the annihilation operators of the 
Bunch-Davies initial state in region II (where m^ = 2), determined in eqn. (6.26) to 
be 



lCl7r/2.(+)_ -[Ck/2.(-) 

(+) ^ ■^c ^ ■^c 



U] 
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where 

We continue these mode functions into region I. For large -u, (e*^"/ sech [-u]) ^ (l/o") • 
(o"/2)~*^. Therefore in region II near the light cone one has 



^^ '^ 47rv^ ■ cosh[r] ■ a Kl) ' ^"'^^ 



Recall that (as shown in appendix A) the matching conditions across the light cone 



are 



,+.C-l . ^ ^ ,2,) . 5ilM . i+.C-i, 

cosh [r J sinh[4j 

a+<-i ■ ^ > 

coshH ' ^g^^^ 



cosh[r] 



a-<-' . ^^ -. 0, 



.-<-'. ^^ (-2.). !^.r=c-., 

cosh [r J smh[^J 



Therefore, for C > only f^ has a nonvanishing continuation into region I, and 
similarly for C < only /i has a n 
C > 0. In region I near the light cone 



'C 
similarly for C < only /i has a nonvanishing continuation into region I. We set 



.5) 






Ve''^ - e-^*^ 27rv/C sinh[^] 



We now change variables to the gauge-invariant potential $ using eqn. (7.24), 
thus introducing a factor of t^ iTGV^(f)/{^iC + 2). We write $ as a quantum- mechanical 
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operator, so that the s-wave component 



$ 



dc 



$ 



(+) 



$ 



(-) 



sin[a] 
sinh[^] 



^.6) 



where 



$ 



+ 



nGV,,- 



.</2 



1 



2^l^/( Ve'^C - e-'^C 



2^C 



-< + 2 



F(t; C) ac + e 



-</2 



2-iC 
< + 2 



F{t;Cyd_ 



17) 



and F(t; (^) satisfies eqn. (7.20) and is normahzed to behave as F{t; C) ~ ^ for 

small t, so that 



F{t; C) = 2(1-*'^) e+'<'^ e^ 



c- 



3(C 



$*^+'' is the 'positive frequency' part of $ — that part which annihilates the Bunch- 
Davies vacuum. The expression for $^ (which consists of one quantum mechanical 
degree of freedom) involves two harmonic oscillator degrees of freedom, rather than 
just one such degree of freedom, because there exist correlations between regions I 
and 11. (If one considers the Bunch-Davies vacuum restricted to region 1, one has a 
mixed state rather than a pure state because of correlations across the light cone.) 

As discussed in the previous section, F{t; () oscillates for small t, while the mode 
labelled by ( is within the Hubble radius, and then for large t, when the mode is well 
outside the Hubble radius, F(t; Q ^ C\^ + C\ e^*. We are primarily interested in 
the growing mode at late times. Therefore we compute the power spectrum for $ as 



lim($^(t)$C,(t)>=P<j>(C)5(C-0 



^$(0 



GVA 4 coth[7rC] 



H 



(8.9) 



3 ae + iy 



where by dimensions we have restored the Hubble constant H during infiation. In 
this formula we maintain ( as a dimensionless variable. In terms of the Laplacian on 
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hyperbolic space, we have — V^ = k'^ = H'^{('^+1). Using x = 167rG'(V^/V|)$ (during 
inflation), one obtains that the power spectrum for x (which has a normalization that 
relates more directly to the density perturbations seen after reheating) is 

With the conventions used here P ~ C "^ corresponds to scale invariance. This is 
seen for example by computing ($^(0)) using the small ^ limit of eqn. (8.6) and 
noting that there is a logarithmic divergence at large (. This result generalizes the 
standard calculations ' ' ' of fluctuations produced in a flat universe during in- 
flation to our open inflationary scenario. We compute the CMB perturbations and 
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other cosmological consequences in a separate letter. 

Our result differs from that of Lyth and Stewart , and of Ratra and Peebles 
who assume different initial conditions for the quantum fields, only by the factor 
coth[7rC/2]. This is very close to unity on all scales accessible to observation (e.g. in 
the large angle CMB anisotropy), so we expect very similar phenomenology to the 
Lyth-Stewart-Ratra-Peebles spectrum. However, it is probably fair to say that our 
calculation, in which the initial conditions are physically justified, puts the result on 
a firmer footing. 



9. Concluding Remarks 

Finally, we conclude with the following comments: 

1. Although it is possible, as we have shown, to extend inflation from a theory 
that predicts fi = 1 to one that predicts Q < 1, we do not believe that a similar 
extension is possible for fi > 1. To be sure, one can construct a closed inflationary 
model by postulating positive spatial curvature at the beginning of inflation, but 
then inflation no longer solves the smoothness problem, and much of the original 
motivation for inflation is lost. 
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The fundamental distinction between closed and open inflation is that negative 
curvature can be introduced locally, by an event localized in space and time (in our 
scenario the nucleation of an isolated bubble) which propagates at a speed asymp- 
totically approaching the speed of light, thus producing smooth surfaces of constant 
negative spatial curvature. By contrast, to single out surfaces of constant positive 
curvature subsequent to the beginning of inflation would seem to require some sort 
of nonlocal process. 

2. One possible source of error in our calculation of the density perturbations 
is the neglect of gravitational perturbations in region II. In region II we used the 
"stiff' approximation, in which the effect of metric perturbations on the scalar field 
perturbations is ignored. We know that the 'stiff' approximation works well while 
modes are well within the Hubble radius. Therefore, for (^ ^ 1, in other words for 
modes well within the Hubble radius at the beginning of the region I phase of inflation, 
we do not expect using the stiff approximation in region II to be a significant source 
of error. However, modes with (^ ^ 1 are never within the Hubble radius at the 
beginning of region I. This is because as one approaches the beginning of region I, 
formally at least, the universe becomes curvature dominated, and the size of the co- 
moving Hubble radius approaches a constant. Fortunately, since consistency with 
observation requires that Q > 0.1-0.2, all scales accessible to observation (i.e., within 
the present apparent horizon) are never far from the range of validity for the stiff 
approximation in region II. 

3. In calculating a density perturbation spectrum, for purposes of computational 
simplicity, we assumed that the effective mass of the infiaton field changes instanta- 
neously from m^ = 2H^ outside the bubble wall to m^ = inside the bubble wall, 
which is assumed to be infinitely thin. One would expect a slightly different spec- 
trum for a wall of finite thickness, especially for large wave numbers. In particular 
for ( large in relation to the inverse of the wall thickness, the modes respond to the 
change in mass adiabatically rather than according to the sudden impulse approxi- 
mation. We have extended our calculation to the case m ^ 2H case, in which case 
the 'thin wall' approximation would be fully justified. The results will be presented 
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elsewhere. On scales accessible to observation it turns out that varying m? / H"^ has 
little effect on the power spectrum. 

4. Refs. 7 and 8 consider a scenario of open inflation restricted to region I. As 
an initial condition they impose the requirement that the initial state is annihilated 
by the operators associated with modes that have positive frequency with respect to 
conformal time (i.e., with the asymptotic behavior t~*l^l~i near the light cone). This 
choice of initial condition is not connected to what happens outside the light cone 
(prior to the coordinate singularity at t = 0) and not surprisingly gives unphysical 
behavior near the light cone: It is a state akin to the Rindler vacuum for ordinary 
Minkowski space. For such a state the stress-energy observed by a freely falling 
observer crossing the light cone diverges; for the Bunch-Davies vacuum there is no 
such divergence. 

Except for the factor of coth[7r,^] in eqn. (8.10), our power spectrum is in agree- 
ment with refs. 7 and 8. This is precisely the expected discrepancy, resulting from 
the Bogolubov transformation, because as far as the evolution in region I is concerned, 
we are in agreement with refs. 7 and 8. 

5. In the model of open inflation scenario presented here the "big bang" singular- 
ity at t = in the open expanding FRW universe is not a genuine singularity. Rather 
it has been reduced to a coordinate singularity, similar in character to the coordinate 
singularity of the black hole horizon in conventional Schwarzschild coordinates. A 
freely falling observer passing from region II into region I would not experience any 
singularity. With this 'big bang' singularity removed, it is tempting to contemplate a 
universe eternal in the backward time direction, with no initial singularity or Planck 
era. Unfortunately, as shown by Vilenkin and by Borde and Vilenkin, eternal in- 
flation backward in time seems to be inconsistent with a finite bubble nucleation 
rate. 

6. We did not consider the generation of gravitational waves in our model. One 
would expect a sizable contribution to the CMB from tensor modes when the energy 
scale of inflation is close to the Planck scale, just as in flat inflation. 
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Note Added 

We have recently extended the calculation of the spectrum of density perturba- 
tions to the case of arbitrary false vacuum mass m (assumed positive). The result 
is surprisingly insensitive to the precise value of m, changing little as m varies from 

Too] 

zero to infinity. We have also received a recent paper by Yamamoto et al.' in which 
they attempt to perform the same calculation as we have done here using instead 
analytic continuation of the Euclidean vacuum modes, with no change in the mass 
of the scalar field occurring across the bubble wall (they assume m^ is constant and 
much less than H^ everywhere in de Sitter space). If we make this assumption (which 
seems hard to justify physically) in our approach, our result is different from theirs. 
We attribute this to the failure of the Euclidean continuation method to correctly de- 
scribe the matching of quantum field modes across the bubble wall (this is discussed 
in detail in ref. 37). 
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APPENDIX A — Matching Conditions for Minkowski Space. 

In this appendix we present the details of the derivation of the matching conditions 
for Minkowski space for the expansions in regions I and II assuming a massless scalar 
fleld. The calculation is carried out only for the s-wave. 

Our strategy is the following. We expand in terms of the more customary mode 
expansion 



oo 



0(rm,tm) = j dk k^JQ{krm) x i^a^'^k) cos[A;t^] + a^''\k) sin[A;t^]| , (A.l) 



which is nonsingular on the light cone and valid in both regions I and II. 
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The expansion in region II 



0(^,r) = f dQs a+'<''-' X {41{Cds)Q^'\T;Cds)+AfhCds)Q^''\r;Cds)} 



— oo 

+00 



r+iCdS-1 



rfCrfS^^^I^ X [a\]\Qs) cosiQsr) + 4"ACd5) MCdsr)} 

—00 

(A.2) 

contains enough information to determine the coefficient functions a^^'{k) and a^"''{k), 

because region II contains a Cauchy surface, defined by r = 0. 

It follows that 



00 



+00 00 / (s) / \ 

—00 \ jj 



MW(A:;Crfs) J' V4/(Crfs) 



where 

CX) 

o r 

-iCds 



M^'\k-Cds) = ^jdasm[ka]a+ 


= ^r(l + ^C.5)A;-^--'coshf^ 



and similarly 



(A.4) 



M('*)(A;;Cds) = ^ f dasm[ka] a+^'^^s-i 

"" (A.5) 
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where we have used the relation 



dt t" sin(at) = a"^"" 1(1 + a) cosh 



iTTO; 



(A.6) 



We now compute the coefficients of the region I expansion 



^^^''^ = j^^' ^^ ^ {^^^^(^^) t+^^-V 4-^(0.) t-^-^ 



(A.7) 



in terms of a*-*-' (A;) and a^"'>{k). 
The relation 
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where t > is arbitrary but fixed, may be rewritten as 
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Taking into account eqns. (A.3)-(A.5), we may write 



AriCh)=t 



=FiCh+l 



-<.Ti) ^4 



X - — / dC sinh^f — 



sinh[^] 



r, sin f A;t sinh ^) f 2 , , -r o 

A;tsinn[^J J vr 

-oo 

X I {iCds)cosh (^\ COS [H cosh e] 4/(Crf5) 



2 



+ (tCds) sinh f^^ sin [H cosh e] 4/ (Cds) 
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We evaluate 
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first carrying out the integration over k 
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and then the integration over ^ 
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one finally obtains 
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Note that 
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vanishes where one of the two delta functions S{(h ± (ds) is nonzero; therefore, eqn. 
(A.14) becomes 
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From the well-known relation r(a;)r(l— x) = tt/ sin(7ra;), it follows that r{ix)r{—ix) 
7r/(xsinh[7ra;]), hence the last line. 
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Similarly, we evaluate 
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Integrating over k, we obtain 



dk k^"^-"'-^ 



2 r(-<ds) cosh 2 



cos{kte ^) — cos(A;te~''^) 
TrCdS 



(te" 



C^+iCdS _ Ug+i\+KdS 



'-i) T{-iCds) cosh 



<rf5\ ^+^c,s siniQsO, 



so that we have 
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Therefore, for C > 0, 
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We may summarize this result in a more compact and more intuitive way as 
follows. In terms of the mode functions, as one passes from region II to region I, one 
has (for C > 0) 

(2.) . ^ . t+^^-\ 
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(A.21^ 



sin[C^] ^„i(^_i 



cosh[r] sinh[,^] 

We finally point out that these matching conditions are generally applicable to curved 
space and nonvanishing mass. It is the behavior of the leading singularities at the 
light cone that determine the matching conditions; therefore, introducing mass and 
spacetime curvature gives subdominant effects. 
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Figure Captions 

Fig. 1 5*0(3, 1) symmetric coordinates. Spacetime is divided into three regions — 
regions I, II, and III. In the zeroth order solution the scalar field is constant on 
each of the hyperboloids, and on the light cone of the origin O. The dashed 
line shows the origin of spherical region I coordinates (.^ = 0). 

Fig. 2 Slice through a bubble. In the (lower) Euclidean domain a bubble nucleates 
via the Coleman-De Luccia instanton. In the (upper) Lorentzian domain the 
bubble expands classically at a speed approaching the speed of light. Region I, 
the bubble interior, is an expanding FRW universe with spatial hypersurfaces of 
constant negative curvature. Region II is the exterior of the expanding bubble. 
The scalar field, and background density are constant along curves shown. The 
future null cone of the center of the nucleating bubble is shown as the line r = t. 
The background scalar field takes the value 0n on this null surface. 

Fig. 3 Potential for 'Old'+'New' Infiation. Initially the infiaton field is stuck in the 
false vacuum F, during the epoch of old infiation. This is exited by the nu- 
cleation of a single bubble, which eventually grows to encompass our entire 
observable universe. Instead of tunneling to the true vacuum T, the scalar field 
tunnels onto a slow rollover potential, at the point (pn, and the interior of the 
bubble expands quasi-exponentially, during the epoch of new infiation. In our 
scenario, this is shortened compared to conventional infiation so that the nega- 
tive spatial curvature survives. Then as the field begins to roll more rapidly, and 
starts to oscillate about the true vacuum value (pr, reheating occurs, converting 
vacuum energy into radiation and matter. 
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